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'~T' Abstract 

O The stability of the fundamental defects of an unstretchable flat sheet is examined. 

• This involves expanding the bending energy to second order in deformations about the 

Cd defect. The modes of deformation occur as eigenstates of a fourth-order linear differential 

^ operator. Unstretchability places a global linear constraint on these modes. Conical 

F^ defects with a surplus angle exhibit an infinite number of states. If this angle is below a 

C critical value, these states possess an n-fold symmetry labeled by an integer, n > 2. A 

^ nonlinear stability analysis shows that the 2-fold ground state is stable, whereas excited 

I— I states possess 2(n — 2) unstable modes which come in even and odd pairs. 

oo 1 Introduction 

o 
o 
in 

[^ In general, when an elastic sheet bends it will also stretch. While the energy associated with 

O bending is purely geometrical-depending on the surface curvature-the energy associated with 

^ stretching will depend on additional material degrees of freedom. However, if the sheet is very 

thin compared to a typical radius of curvature, stretching will typically become far more costly 

than bending; thus, if stretching does occur, it will tend to be confined within a set of sharp 

rN peaks connected by a network of folds in the sheet (Tj [21 E]. This observation has attracted 

c3 many physicists in the last decade, who have started to look at the properties of these defects 

in detail (a nice review is provided in reference [4J). The limiting behavior of thin sheets has one 

striking feature: it involves only the geometrical degrees of freedom of the sheet, in particular 

the bending energy depends only on its shape; unstretchability constrains what shapes are 

accessible. This simple picture turns out to provide a description of the behavior of a range of 

real materials on widely differing scales: it is equally valid for a sheet of graphene on microscopic 

scales as it is for a sheet of paper on macroscopic ones [H |5]. Corresponding materials of soft 

matter have been studied only recently, mainly inspired by biology. It was shown for example, 
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that growing tissues or thermo-responsive gels can be described within the same limit if they 
are sufficiently thin [51 [71 |H]. A general treatment of their stability is, however, still lacking. 
The goal of this paper is to analyze the stability of the basic building block of such materials, 
the conical defect. 

Approaching the stability of the equilibrium states of an unstretchable sheet is complicated by 
the local nature of the constraints. These constraints distinguish their treatment from that 
of a fluid membrane which can also be modeled geometrically [9J. In fluid membranes the 
penalty associated with shear is very small so that it behaves like a fluid along tangential 
directions. If there are constraints in the latter they tend to be global in nature: area, enclosed 
volume or-in a bilayer-area difference between inner and outer layer may be constrained [lOj . 
In an unstretchable sheet, however, the constraints are local. Unstretchability translates into a 
constraint on the metric tensor induced on the surface [TT]. While the two limits lend themselves 
to a geometrical description the physical behavior they display is qualitatively very different. 
This is reflected in how one goes about examining this behavior. 

The metric constraint lends rigidity to the surface with respect to its behavior as a fluid mem- 
brane. Whereas a shape that is stable as a fluid membrane will also be stable if it is un- 
stretchable, the converse is not true: unstable modes of a fluid membrane will tend to become 
inaccessible when the constraint on the surface is accommodated. As a case in point, there 
are conical defects that are generic stable features of unstretchable flat sheets. If the apex 
of the cone is part of the surface it is not possible to smooth it out by surface deformations 
that preserve the conical metric. In a fluid membrane, however, the energy will be lowered by 
smoothing out the singularity at the apex. Thus, conical defects will never appear as stable 
features of a fluid membrane. In this paper the discrepancy will be quantified. 

The stability of simple fluid membrane equilibrium states-spheres and cylinders-was addressed 
by Helfrich and Ou-Yang in the late 80s [T2]- This involves the evaluation of the second variation 
of the bending energy about these surface states. This was a non-trivial calculation and, not 
surprisingly, their approach was adapted specifically to these states. Later a framework was 
developed which does not depend explicitly on the equilibrium state [13j. 

To accommodate the unstretchability of the sheet, we will adopt the method of Lagrange 
multipliers in the calculus of variations to enforce the corresponding metric constraint. Like the 
construction for fluid membranes developed in [13J the framework we will develop will not depend 
explicitly on the particular equilibrium shape. Our focus will be on flat sheets, which can be 
modeled as surfaces with vanishing Gaussian curvature, and we will derive an explicit expression 
for the second variation of the bending energy subject to this constraint. The unstretchability 
manifests itself in the important respect that the second variation, unlike its fluid membrane 
counterpart, depends now on both the normal and tangential deformations of the surface. This 
is, as we will see, because the constraint spoils the identification of tangential deformations 
with reparametrizations. That an exact framework exists at all is somewhat surprising; that the 
calculation is largely tractable analytically is even more so. 



This geometrical framework will be used to examine the stability of the equilibrium states 
of a circular flat sheet with a conical defect at its center. Such surfaces are in one-to-one 
correspondence with closed curves of fixed length on the unit sphere. We will focus, in particular, 
on a sheet with a surplus angle at this point |Zl[B] (see also [14] for a related example involving 
nematic solids). Unlike their counterparts with an angle deficit such cones possess an infinite 
number of equilibrium states in the absence of external forces. If the surplus angle is sufficiently 
small, so that self-contacts do not occur, these shapes are characterized by an n-fold symmetry 
beginning with a 2-fold. To assess their stability we obtain an explicit expression for the second 
variation of the energy about each of these states. This will involve a self-adjoint fourth-order 
linear differential operator, £, which depends explicitly on the curvature of the conical state. 
The accessible modes of deformation are described by the eigenmodesof this operator consistent 
with a global linear constraint associated with the unstretchability. A negative eigenvalue signals 
an instability. If the surplus angle is small this constraint is tractable analytically and an exact 
Fourier analysis of stability is possible. 

The rotational invariance of the energy implies the existence of three zero modes of deformation 
for each equilibrium state of the conical defect. An exact expression for these modes will be 
provided. Two of them correspond to rotations about axes orthogonal to the axis of symme- 
try. For small surplus angles they possess a pair of nodes. The third zero mode is odd and 
corresponds to a rotation about this axis; it has 2n nodes in this limit. All eigenmodes with an 
intermediate number of nodes will be unstable. 

In general, an exact treatment of stability eludes us. So we develop an approximation for larger 
surplus angle by introducing a Fourier expansion of the modes; using the Gram-Schmidt process 
we construct an orthonormal basis satisfying the global constraint. The matrix elements of the 
linear operator are evaluated. While the ground state is stable-like its ice-cream counterpart-as 
one would have predicted, all excited states are unstable. The detailed nature of this instability 
is, however, rather interesting. The n-fold beginning with n = 3 will possess 2(n — 2) unstable 
modes. These modes occur in pairs with even and odd parity with respect to the n-fold. For 
low surplus angles, each mode consists of a single frequency. As the surplus angle increases 
higher frequencies enter and the number of nodes of each mode increases. It is nevertheless 
possible to order modes in a manner analogous to the small surplus angle limit. The spectrum 
of the operator C will be discussed and the dominant mode of instability will be identified for 
each excited state. 



2 Second order variations with a metric constraint 

We will model the unstretchable sheet in terms of a surface which is described parametrically by 
three functions X(u^,u^), X = (X^, X^, X^) providing its position in three-dimensional space. 
The bending energy associated with the folded state of the surface is quadratic in curvature, 

Hb[X] = ]^J dAK\ (1) 

3 



where dA is the area element induced on the surface and K = Ci + C2, where Ci and C2 are 
the two principal curvatures. Hb[X] thus involves only the geometry of the surface. Once the 
geometry is specified, the bending energy is also. We will work in units in which the rigidity 
modulus is unity. 

The induced metric tensor on the surface is given by gab = e^ ■ e^, where e^ = daX, a = 1,2, 
are the two tangent vectors of the surface adapted to the parametrization by u^ and u^. This 
parametrization will be fixed once and for all. Geometrically, the unstretchability of the sheet is 
then the statement that the surface is isometric, i.e., the only admissible deformations are those 
keeping gab fixed. A natural way to keep track of this constraint in the calculus of variations is 
to replace Hb[X] by the functional [llj 

H[X, T^"] =Hb[X]-\J dAT^' {gab - g^^) • (2) 

The symmetric tensor T"^ is the set of Lagrange multipliers associated with the constraint that 
gab coincides with some fixed metric g\^^ . In particular, we will be interested in a surface that 
is a circular flat disc with a conical singularity, described by the line element 

df = dr^ + r^ds" , (3) 



where r is the distance to the apex of this cone and s ranges between and some value 27r + (y9 
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However, our framework-at least in its initial stage-will not require us to specify g)^/ explicitly. 



Recall that a surface is completely determined once the metric tensor gab and the extrinsic 
curvature tensor Kab = e^ ■ dbU are specified [T5l US]. The symmetric second rank tensor Kab 
is a measure of how fast the normal vector n rotates into one direction as it is moved along 
another. Unlike gab, it depends explicitly on this vector. The principal curvatures, Ci and C2, 
are the maximum and minimum values of the shape operator Kab- They are assumed along 
two tangents-the principal directions-which are always perpendicular. Fixing the metric places 
a constraint on deformations of the surface that is more severe than may, at first, appear. A 
familiar example is provided by the pingpong ball in which curvature and spherical topology 
together conspire to produce a geometry that is rigid under isometric deformations. 

In order to say something about stability we first identify the Euler-Lagrange equations that 
describe the equilibrium states subject to the constraint. Consider a deformation X — > X+5X, 
rpab _^ T"-'' -\- ST"-'' , As described in reference [TT], where the variational framework developed in 
references [17j and [18J was adapted to accommodate the metric constraint, the corresponding 
change in H at first order is given, modulo a boundary term, by 

6H= I dA (S^n + S^^eb) ■ 5X - ^ 6T"\gab - g^^) , (4) 

where 

(5) 

(6) 



J 


/ dA {£^n + S^eb) ■ 5X - 


- - 5T"\gab - g, 


s± 


= -V^K+^K{AKg- 


-K^)~KabT"'' 


^1? 


= VbT''^ 





Here Kg = C1C2 is the Gaussian curvature. The constraint adds a term linear in curvature to 
the normal Euler-Lagrange derivative £±. Unlike a fluid membrane, the tangential counterpart 
is not identically zero. The equilibrium of the sheet is then described by the Euler-Lagrange 
equations: 

8^=0, 8^ = 0, gab = 9^S. (7) 

In equilibrium, the tangential stress T"^ associated with the metric constraint is conserved. 
This equation involves only the intrinsic geometry. The extrinsic geometry enters through the 
boundary conditions on the free surface of the sheet. In general, these conditions will involve 
the introduction of a boundary layer in which stretching is admitted fH E]. Fortunately, in 
the case we will treat here, the energy itself will be independent of the precise nature of these 
conditions. The only component of T*^^ which enters the stability analysis is the trace T^-^Kah 
coupling to curvature in 8^. In [llj it was shown that this component is insensitive to boundary 
conditions. 

We now need to extend this framework to address the stability of equilibrium shapes. To do this 
it is necessary to expand H to second order about solutions of the Euler-Lagrange equations. 
Note that when the Euler-Lagrange equations are satisfied. 



5'H = dA 58^ (n ■ 5X) + 58^^ (e^ ■ 5X) . (8) 

Here 

88^ = 6x8^-Kai,ST^'; (9) 

68^ = Vb5T''\ (10) 

where 5x8i_ represents the first order change in 8a_ given by Eq.([5]) under the isometry 5X, 
maintaining T"-'^ fixed. Let us examine 58^-. Note first that terms involving the deformed 
Christoffel connection F^^ vanish whenever 5gab = 0. A more general statement is possible: 
whenever 5gab = 0, the variation commutes with covariant differentiation. We now integrate 
by parts to express the second term appearing in Eq.(|8|) in the form 



Vb^T'^" (e, ■ 5X) = ST'^'Kab (n ■ 5X) - ^T"" (e„ ■ V,<5X) , (11) 



modulo a divergence. The first term cancels an identical term appearing in 68±_ (see Eq.([9 
the second term vanishes because 

5gab = ea ■ VbSX + e, ■ VaSX = . (12) 

What is left is the surprisingly simple expression 

S'^H= [ dASx8±{n-SX). (13) 



Superficially, it might appear that the second order deformation is identical to the one that 
would have been obtained (see reference jl3j) had T"^ been treated as a parameter. All that 



one seems to have achieved is the satisfaction of having dotted all of the i's, confirming that 
what one might have guessed is also correct. The presence of the constraint, however, spoils the 
identification of tangential deformations with reparametrizations of the surface: in particular, 
a tangential deformation of £± does not vanish. While £± is a surface scalar, it involves 
the tangential stress T"^ which does not depend explicitly on X. It is thus not completely 
determined by the local surface geometryj^ As a consequence 6x£± is not the same as 6±S± 
in an equilibrium state, and it is no longer legitimate to replace one by the other as it was in 
the case of a fluid membrane (see, for example, [13]). 

While this does not appear to bode well, the evaluation of 5x£± turns out to be simpler 
than that of 5±£±. This is because it is possible to exploit the fact that isometric variations 
commute with covariant differentiation. At a later stage in the calculation, it will be necessary 
to implement the metric constraint explicitly to eliminate tangential deformations in favor of 
normal ones. 



To obtain Eq.(13) we discarded a boundary term in Eq.(ll) involving ST"-''. Had we kept track 
of it, we would have seen that it cancels a counterpart originating in the boundary term that 
tags along on the right hand side of Eq.(|4]) when derivatives of SX are collected in a divergence. 



There remains to determine 6x£± appearing in Eq.(13) 



It is useful to think of the normal Euler-Lagrange derivative, £±, defined by Eq.([5]) as a sum 
of two terms. The bending part depends on the extrinsic geometry only through the mean 
curvature. Using the fact that both the Laplacian and the Gaussian curvature are invariants 
under isometry, it is simple to see that 



6x£± = -W'5K + 2Kg SK - -K'SK - T^'SKab ■ (14) 

The bending contribution only requires the evaluation of 5K. The tensor SKab does, however, 
enter the source term associated with the constraint. In appendix |A| 5Kab and 5K are given 
in terms of the deformation 6X. In the following we will specialize the obtained expressions to 
study the stability of conical defects on a flat unstretchable sheet. 

3 Conical Deformations 

3.1 Cones as curves on spheres 

We will use the representation of the conical defect as a closed curve F : s i— )■ u(s) on a unit 
sphere (see Fig. [I]). If r is the distance to the apex, then a parametric description of the cone 
is provided by the identification X(r, s) = ru{s). It is convenient to parametrize the curve 



^A small tangential deformation of any geometrical tensor whose spatial dependence is determined locally 
by X can be identified with a reparametrization of the surface (a Lie derivative). It thus vanishes when its 
argument vanishes. 
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Figure 1: The geometry of the conical defect. 

by arc-length. The total length of the curve Se = 27r + ip^ will be invariant under isometry 
Indeed, the intrinsic geometry of the cone is characterized completely by s^. For technical 
simplicity let us suppose that this cone has a fixed finite radius R. If s^ = 2n, this cone will 
be isometric to a planar disc. If Se < 27r, it will exhibit an angle deficit at its apex. The 
unique equilibrium configuration in the absence of external forces is then an axially-symmetric 
ice cream cone (see again Fig. [I]). Our focus will be on cones with a surplus angle s^ > 2n. 
Unlike their counterparts with an angle deficit, as described in [8J, there will be an infinite 
number of non-trivial equilibrium states. In particular, if the surplus angle is sufficiently small 
so that self-contacts do not occur, these states will be symmetrical with an n-fold symmetry 
labelled by an integer n > 2. The axis of symmetry will be aligned with the basis vector k of 
the Euclidean coordinate system {i,j, k). Once self-contact occurs the symmetry of the cone 
is broken and skewed geometries emerge [19J. 

If t is the tangent vector to the curve w(s), and n = uxt, then {u,t, n} forms an orthonormal 
frame satisfying 

u = t , t' = —K,n — u , n' = K,t , (15) 

where k is the geodesic curvature of the curve on the sphere (which, confusingly, turns out to 
be its normal curvature if it is viewed as a curve on the cone). 



The extrinsic curvature tensor is given by 



K„h = r 




K 



(16) 



The flat direction is along u. The Gaussian curvature Kg vanishes. The trace K = n/r 
diverges at the apex of the cone. This translates into a logarithmic divergence in the bending 
energy at this point. 



The normal projection of the Euler-Lagrange derivative is given by 



£i 



1 



K 



-K^ + (1-C||)k 



:i7) 



where C\\ = —r'^T^^ is a constant related to the stress along the tangent direction t [llj. The 



'\ 
vanishing of £^_l determines the equilibrium shapes of the cone. Its first integral 

1 



r 



f<2 



4 "^ 



[li 



can be solved in terms of elliptic functions where J = Jk \s a constant vector aligned with k 
which is associated with the rotational invariance of the conical bending energy [llj. It is given 
by 

~2 



Ci 



n 



K,'t + KU . 



(19) 



Solving the first integral (18) for a cone with surplus angle ipe and n folds one obtains [8] 



k{s) 



A^^K [k] 
'S 



sn 



2slC[k] 



(20) 



where S := (27r + Lp^)/2n = Se/2n. Closure of the surface in Euclidean space determines 



the parameter k (see the discussion following Eq. (44) in section 4.2). It is directly related 

1 -4(1 + k)IC[kf /S^ and 

Cn — GAklC [k] /S^. In Eq.(20) the function sn (s, k) is the sine of the Jacobi amplitude 



to the stress Cy and the constant of integration J via C| 
J= ^~ 



am {s, k) with parameter k. The symFol /C [k] denotes the complete elliptic integral of the first 
kind t2Qj. 



3.2 Deformed cones 

To assess the stability of the cones with surplus angle, we need to examine the second variation 
of the energy. Unlike the unique ice-cream cone, the question of stability is non-trivial. For 
technical simplicity we will limit the discussion to cones which do not exhibit self-contact, i.e. 
Lfe < 7.08 for n = 2, for example [8j. 

The local deformations of the cone will not only be constrained by isometry they must also 
preserve its apex. The deformed cone is thus another cone. This can be described in terms of 
the deformation of the trajectory on the unit sphere, u{s) — )■ u{s) + 5u{s). The only non- 
vanishing component of the deformed extrinsic curvature tensor is then given by SxKss = rSuH. 
One thus only needs to evaluate 5uK. 

The unstretchability of the cone translates into the statement of conservation of arc-length on 
this sphere. In particular, the angle surplus is preserved. The change in arc-length induced by 
5u is given by 6ds = t ■ {5u)' ds. Thus fixed arc-length is captured by the identity 



t ■ {6uy = , 



(21) 



It is useful to recast Eq.(21) as a relationship connecting the two non-vanishing projections of 
5u, ^ = n ■ 6u and ^ = t ■ 5u: 

^' = -K$. (22) 
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This is a special case of Eq.(56) with the identification = r$, ipa = r^t ■ Ba- In particular, 
note that unless k vanishes, every non-vanishing normal deformation will be accompanied by 
a compensating tangential deformation. In general, a global constraint on $ is implied by 



Eq.(22), namely 



ds K^ 



(fi:|$) = 0. 



(23) 



Because arc-length is conserved, derivatives with respect to s commute with the variation. 
Thus, in particular, 

6t = {6uy . (24) 

Sn is thus expressed in terms of Su as follows 

Sn = —{n-5u)u—{n-5t)t 

= -{n-6u)u-{n-6u')t. (25) 

We are also now in a position to write down how the curvature k = —n ■ t' changes. We have, 
for a length preserving deformation, 

^ui^ = —fi ■ [St)' — 6n ■ t' 

= -n ■ [6u + 6u"] . (26) 



The second structure equation appearing in (15) was used on the second line. Casting 6k in 
terms of the projections, $ and ^, we obtain 



Sk = -$" - $ + [k^]' . 



(27) 



This equation is a special case of Eq.(58) 



3.3 Second Variation 



The second variation (13) of the energy now reads, modulo a constant multiplicative factor 



originating in the integration with respect to r. 



6^H = - lds^\5K" + (-k^ + I-CASk] 



(28) 



This is also what one would have guessed without going through a detailed analysis. The naive 
term linear in the variation of the multiplier 5C\\ (assumed to be constant) vanishes on account 
of Eq.(23). The justification for ignoring 5C\\ is found in section |2| Modulo this caveat, the 
calculation adapted directly to the conical geometry does get the answer right. 



Let us now substitute the expression for 5k given by Eq.(27) into Eq.(28). One finds that the 



dependence on the tangential deformation \E' can be collected in a total derivative so that the 
second variation can be expressed completely in terms of <l>: 



5^H= i ds<^i<^"" + Vi{s)<^" + V2{s)<^\ 



(29) 



where 



Vii-s) = ^K^ + 2-Cii 

= l^' + [^ - (1 - cii)]'^' + liJ' - ci) + 1 - qi • (30) 



The second presentation of V2 uses the Euler-Lagrange equation as well as the quadrature (18). 
The details of this derivation are collected in appendix [B] 

The obtained expression is still not in a manifestly self-adjoint form. Rearranging derivatives 

$(^1$')' = Vi<!><!>" + ^ F/(<l>2)' , (31) 

permits Eq.(29) to be cast in the form 

S^H = (ids<^[^"" + [Vi{s)^']' + fv2is) + lvl'{s)\ $} . (32) 



We write 



H= ids^C^= ($|/:|<l>), (33) 

where the fourth order linear differential operator C, given by 

is self-adjoint. As such, C is guaranteed to possess real eigenvalues. If all these eigenvalues are 
positive (for a given n and Lp^), the corresponding equilibrium shape will be stable. 

4 Stability 

4.1 Small Surpluses 

Before analyzing the spectrum of C in the full non-linear theory, we provide a complete stability 
analysis of cones with a small surplus angle (/?e- In this case k is also small, so that the 



quadrature (18) implies that the equilibrium curvature is given hy k = Kosin(ras) where n > 



2 G N (we fix the phase without loss of generality) which identifies Cy = 1 — n"^ and J^ = Cj 
to lowest order. The ground state corresponds to n = 2. In this approximation, both Vi and 
V2 are constants, Vi = n'^ + I and V2 = n^ so that 

C=— + il + n^) — +n\ (35) 
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Figure 2: Eigenvalues Am in the small surplus limit as a function of mode number m for 
n = 2 (circles), n = 3 (squares), and n = 4 (diamonds). The joining lines correspond to 



the eigenvalue Eq.(36) but have no physical significance for non-integer m; they have been 



added as a guide to the eye. Note that for m > 1 each point corresponds to two modes 
due to the two-fold degeneracy of the small surplus limit. 



Periodicity implies that the modes of deformation are represented by a constant (m = 0), 
sin(ms) and cos(ms), m > 1 G N. For fixed n, the eigenvalues of C are then labeled by the 
integer m = 0, 1, 2, . . . , given by 

X^ = {m^ - r?){m^ - I) . (36) 

In Fig. [2] Am is plotted as a function of m for the lowest values of n. The constant mode with 
m = has positive Aq = n^ ■ All eigenvalues with m>\ possess a two-fold degeneracy. As we 
will see this degeneracy gets lifted for some modes when the angle surplus is increased. 

Four zero modes with A^ = occur at m = 1 and m = n. The mode cos{ns) ex k' describes 
a small rotation of the cone about the axis of ra-fold symmetry k. The two modes with m = 1 
describe rotations about an orthogonal axis. Thus these three modes are the zero modes 
predicted by the rotational invariance of the conical bending energy. The fourth zero mode 
sin(ns) is inconsistent with the isometry constraint (23). It is also the only conical mode of 



deformation inconsistent with isometry in the small surplus limit. If isometry is enforced, the 



mode does not persist as an eigenmode. As we will see in section 4.4 if the isometric constraint 
is relaxed it becomes unstable for higher surplus angles. 

The 2-fold ground state with n = 2 \s stable. All excited states of the cone {n > 3) are 
unstable for small surpluses. This is contrary to what was claimed in Ref. [SJ. In total 2(n — 2) 
unstable modes of deformation exist corresponding to m = 2, . . . , n — 1. They lie between the 
zero modes at m = 1 and m = n (see Fig. [2]). All modes of deformations with m > n make a 
positive contribution to the energy. 
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The dominant (fastest) modes of instability will be the two modes corresponding to the lowest 
negative eigenvalue. The corresponding value of m is the integer closest to mo = ^/{l + 'n?)/2. 
For n = 3, the only possibility is m = 2; for n = A, the m = 3 mode has the lower eigenvalue. 
For higher values of n, one would expect a cascade of instabilities. For example, the instability 
associated with the state with n = 7 will proceed by 7 — )■ 5 — )■ 4 — )• 3 — )■ 2, at least for small 
ife- However, predicting the dynamics of this cascade is beyond the scope of this analysis. The 
final part of this work will focus on the extent to which the results of this section hold outside 
the small surplus limit. 

4.2 Zero modes in the non-linear regime 

One can use the symmetry of the system to identify the zero modes of the operator C satisfying 
the linearized Euler-Lagrange equation 

£$ = 0. (37) 

As mentioned above, three zero modes are anticipated by the invariance of the conical bending 
energy with respect to rotations about its apex. Consider the infinitesimal rotation 6u = uxb. 
Its normal projection, given by $6 = bt, then satisfies Eq.(37) for each choice of 6J^ While the 



quadratic constraint t^ = 1 connects the three components, this is not relevant at the linear 
level. 

It is simple to check that these three modes are, in fact, consistent with global isometry: 

{K\^h)=h- idsKt^b- idsn' = Q. (39) 



In our analysis of small Lp^ a fourth zero mode proportional to k was identified. However, this 
mode of deformation was seen to be inconsistent with the global isometry constraint and can 
thus be discarded here. 

The three zero modes can be expressed as projections oft onto the Euclidean coordinate system 



(i,j,k). Using the fact that k is aligned with J we obtain with the help of Eq.(19) 



$^ = fc . t = k'/J . (40) 

which identifies the first zero mode with k' up to a normalization (see Fig. [Sjfor a numerical 
example with n = 3). However, we also know that k' satisfies the linearized Euler-Lagrange 



equation (61) 



Sk" +{-k^ + 1-Ci\)Sk = 0. (41) 



^Note that $& is not a solution of the non-self-adjoint equation 

$"" + Vi(s)$"-l-T^2(s)$ = 0. (38) 



While the total energy does vanish, the integrand appearing in Eq.(29l does not vanish point wise. 
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Figure 3: The zero mode $fc for n = 3 and (pc 
(dotted line). 



^ (solid line), vr (dashed line), and 27r 
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(a) 



(b) 



Figure 4: The zero modes $i and $j for ra = 3 and ipe 
and 27r (dotted line). 
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(sohd line), vr (dashed line). 



This is a direct consequence of the reparametrization invariance of the Euler-Lagrange equation 
and is due to the fact that a rotation about the axis of symmetry can be achieved by a 
reparametrization. 

The second and third zero modes are more complicated and their construction requires a little 
more work. We start by expressing the vector u = cosy^sin^i + sirup sin i)j + cos'dk in 
terms of the polar and azimuthal angles on the sphere, p and ^ respectively. Both angles are 



functions of arc-length s. With the help of the first structure equation appearing in (15) it is 
possible to write 

^i = i ■ u' = —Zi sin p + Z2 cos p and $j = j ■ u' = zi cos p + Z2 sin p (42) 

with Zi = (yp'sint? and Z2 = 'i9'cos'i9. From the projection onto k, $/,, one identifies i!)' = 
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—n' /{J siri'd). Using the identity CDS'*? = u ■ k ^^ k/J, the ^ dependence in the zero modes 
can be eliminated to obtain 



Zi = LP 



1 - —r and Z2 



KK, 



jVj' 



(43) 



To replace the cp dependence, consider the projection n ■ k = {u x u') ■ k = 2:1 sin ^ 
(/tV2 - Cii)/J. From this it follows that cp' = {k^ jl - C\\)I\J{\ - k^ I J^)\ = (J/2){[(j2 - 
2C||)/(J^ — K^)] — 1}. An integration with respect to the arc-length completes the elimination 



of the angular dependences in Eq.(|42]). One obtains 

'l6{-k)}C[kf 



ip{s) 



J' - 2Ci\ S 
2J 2JC[k] 



n 



J252 



, am ( ,k I ,fc 



J 



(44) 



where IT is the elliptic integral of the third kind and am(s,fc)is the Jacobi amplitude with 



parameter k which appears in the expression (20) for the curvature k{s). The closure of the 
surface in M.^ implies that ^(se) = 1-n. From this condition the parameter k can be determined. 

In Fig. |4] numerical examples of $j and $j are shown for n = 3. For small ip)^ the modes are 
proportional to sins and coss as anticipated in the small surplus analysis of the previous section. 
This behavior changes significantly as soon as the non-linearity kicks in: higher frequencies blend 
in and further nodes appear which complicates the identification of zero modes from their 
functional dependence. Note that this problem does not arise with $fe. Being proportional to 
k! , it always possesses 2n nodes no matter what value of p)^ is given (see again Fig. [S]). 



The nodal structure of the zero modes suggests one of the difficulties that will be encountered in 
the full perturbation analysis: a simple identification of the number of nodes with an integer m 
in the manner of the small surplus limit will become unfeasible. How one sidesteps this obstacle 



will be discussed in the Fourier analysis in section [4^ Before turning to the full problem let us 
have another look at the constant mode. 



4.3 Constant normal deformations 

For small p)^, we saw that the constant normal deformation represents a stable eigenmode of 



deformation with m = and Aq 



n 



The constant mode will not generally be an eigenmode. 



However, it is possible to show that b H is positive for such a deformation. 



Rewrite Eq.(32) in the form 



b''H= I ds\{^y -V^{s){^'f 



y^^s) + U/^{s) ) e 



(45) 
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The second term appearing in Eq.(45) is manifestly negative. If $ is constant, then [^ 

5^H = ^^ idsV2{s). (46) 

Modulo the Euler-Lagrange equation, 

V2{s) = li^r + 1^' + (i + ^qi) '^^ + (1 - qi) , (47) 

which is manifestly positive for each n for values of C\\ in the range 1 — n^ < Cy < 1. Outside 
this range V2 can assume negative values. However, we have confirmed numerically for the 
lowest n-folds that the integral itself is indeed positive for all (/9e lying outside the regime where 
self-contacts occur; a transition from one n-fold to another cannot be induced by a constant 
normal deformation. 

Now let us look more generally at the spectrum of C in the full non-linear theory. Let £$j = 
Aj$j, i = 1,2,3, .. . ordered by the magnitude of Aj with periodic boundary conditions on $j 



and C given by Eq. (34). In physics one is perhaps more familiar with second order differential 
operators of the form —dg+V{s). In this case the number of nodes in the eigenmodes increases 
with energy. We have seen that the behavior of the fourth-order operator C is very different. 
While the fourth derivative term will always dominate if the mode oscillates sufficiently, as we 
have seen, for modes commensurate with the n-fold, the negative second derivative term can 
lower the mode energy. It is this term which conspires to produce zero modes exhibiting 2 and 2n 
nodes if the surplus angle is small. For all n > 2, it is not unreasonable to expect that unstable 
modes exist in the nonlinear theory as well corresponding to those with 4, 6, ... , 2{n — 1) nodes 



in the small surplus limit (see section 4.1) 



4.4 Fourier decomposition of modes 

The principal technical obstacle to analyzing the stability of a cone is the non-local nature of 
the isometry constraint: the eigenmodes are orthogonal to k as vectors in a Hilbert space. 
The constraint makes it unlikely that one is going to make further progress without recourse 
to some form of approximation. Our approach to the problem will make use of techniques of 
perturbation theory in quantum mechanics. The language and notation we adopt reflects this 
point of view. 

We first define an orthonormal basis of periodic functions {(pi}, i = 1, . . . ,imax, that is orthog- 
onal to K. These will be constructed out of Fourier modes using the Gram-Schmidt process. 
We begin with the Fourier modes 0i = 1, (f)2m = sin fm^j, and (j)2m+i = cos [m^j for 
m = 1, . . . ,mmax. where s^ = 27r + Lpe as before. We normalize them with respect to the scalar 



■^Such a deformation is consistent with Eq.(23). The identity § dsn = follows from the symmetry of the 
quartic potential. 
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product defined on the Hilbert space, 



2-(5al 



Certain Fourier modes will not, however. 



be consistent with the constraint (fi;|0j) = 0. To see this, let us expand expression (20) for 
k(s) into a Fourier series: 



k{s) 



IGnn 



E(-i)'i 



I Q' 



i+h 



sm 



1=0 



where 



q := exp 



7TlC[l/k] 



{2l + l)n 



2-KS 



(48) 



(49) 



Thus, the sine functions with m = (2/ + l)n and / G {0,1,...} are inconsistent with the 
isometry constraint: 



{l^ I 02(2i+l)n) 



ds K, sin 



(2/ + l)n 



27TS 



i7in{-l) 






72i + l 



47m(-l)' 



cosh 



7rK:[i/fc](;+|) 



7^0. 



(50) 



We immediately note a distinction between even and odd modes of deformation of the cone. 
To construct an orthonormal basis consistent with the constraint we use the Gram-Schmidt 
process on the offending modes, defining the new states for i = 2(2/ + l)n: 



(^ I (pi) 
{k, I k) 



-E 



.) 



t'<i 



') 



(51) 



and normalizing them with respect to the scalar product on the Hilbert space, 0j := (pi/ {(pi \ (pi). 
For i 7^ 2(2/ + l)n we can take the normalized Fourier modes (pi = cpi as defined above. 

For the stability analysis we have to determine the eigenvalues and -vectors of the operator C 
for pairs of (n, Lp^). In the following we consider in particular (/^e = gfg , f , vr, and 27r for n < A. 
An approximation for the corresponding eigenvalues and -vectors can be found by diagonalizing 
the symmetric matrix Cij = {(f)i\C\(pj) numerically for a fixed mmax. hj ^ {1, • • • ,2mmax + l} J 
At least the dominant modes should be included to obtain quantitative results, even though 
low values of mmax should already be sufficient to make statements about the stability. This 
is due to the fact that higher modes will always make a positive contribution to the energy on 
account of the leading fourth derivative term in C 



''Certain simplifications can be exploited in the numerical calculation: Note that ((/'cvcn|'C|(?!)odd) = on 
account of the even parity of £. To simplify the calculation of the integrals in Cij , the non-derivative term 



appearing in ( 34 ) is implemented as 



V2 + I:v;' = -U' 



[2-6(i-q, 



3(j2 -cf) + i- qi 



(52) 
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(a) n = : 


I 




m 


V?c ^0 


7r 
360 


TT 

2 


TT 


2n 





4 


4.0 


3.16 


2.43 


1.41 


1 

















2 

















3 


40 


39.64 


9.75 


3.60 


1.16 


4 


180 


178.75 


58.85 59.46 


23.97 1 24.68 


6.08 6.38 


5 


504 


500.80 


182.27 


80.45 


22.74 


6 


1120 


1113.2 


422.8 422.95 


192.46 1 192.63 


56.67 56.73 



(b) n = 3 



m 


V^c -> 





n 
360 


2 


71 


2n 





9 


8.99 


7.47 


6.04 


3.81 


1 

















2 


-15 


-15.04 


-19.10 


-13.98 


-4.90 


3 

















4 


105 


103.80 


19.82 


8.16 


3.28 


5 


384 


380.88 


105.52 


39.18 


10.25 


6 


945 


938.32 


302.74 307.78 


120.15 127.15 


28.56 33.44 






(c) n = 4 






m 


V:?c ^0 


IT 

360 


7r 
2 


71 


27r 









16 


15.99 


13.53 


11.09 


7.12 






1 





















2 


-36 


-36.03 


-37.75 


-30.91 


-13.87 






3 


-56 


-56.30 


-85.63 


-59.75 


-19.73 






4 





















5 


216 


213.01 


34.10 


14.80 


6.4 






6 


700 


693.5 


164.17 


61.57 


19.99 





Table 1: Eigenvalues corresponding to the modes m of the flat limit for increasing values 
of ife and n = 2, 3, and 4. Each m corresponds to two eigenmodes except m = (the 
constant mode) and m = n (the third zero mode). The two-fold degeneracy persists except 
in those cases where two values are given. 
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To be on the safe side we have taken into account the first two offending modes for (^^ = -^^ 
(by setting mmax = 9 (for n = 2), 14 {n = 3), and 19 {n = 4)), whereas one additional 
offending mode was included for all higher values of ^p^ (by setting m^ax = 13 (for n = 2), 20 
(n = 3), and 27 (n = 4)). 




-0.25 
-0.50 




-0.25 
-0.50 



(a) m = 2, odd mode 



(b) m = 2, even mode 






(c) Unstable equilibrium (d) State n — 3 deformed 
state n = 3 with the odd mode m = 2 



(e) Stable equilibrium state 
n^2 



Figure 5: (a) and (b): odd and even unstable modes for n = 3 corresponding to m = 2 
of the flat limit for V^e = 3I0 (solid line), vr (dashed line), and 27r (dotted line), (c)-(e): 
transition from the unstable state n = 3 to the stable state n = 2 for ip^ = t^. In (d) the 
transparent cone is the initial state n = 3 and the cone with the coordinate lines represents 
its deformation with the odd mode m = 2. 



Tab. [I] presents the lowest eigenvalues for various pairs of {n,ipe). Strikingly, negative eigen- 
values are found for n > 3. A closer inspection reveals in fact that the stability analysis of the 



small surplus limit from section [43] largely carries over into the non-linear regime: whereas the 
2-fold ground state is stable, all excited states of the cone with n > 3 are unstable against 
2{n — 2) modes of deformation (see Fig. [5] for the 3-fold). The mode which was dominant in 
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-0.25 
-0.50 



(a) TO = 2, odd mode 




-0.25 
-0.50 




-0.25 
-0.50 



(b) m — 2, even mode 




-0.25 
-0.50 



(c) m = 3, odd m.ode 



(d) m = 3, even mode 



Figure 6: Unstable modes for n = 4 corresponding to (a), (h) m = 2 and (c), (d) m 
the flat limit for (pe = ^ (solid line), vr (dashed line), and 27t (dotted line). 



3 of 



the flat limit continues to be dominant for higher y^e- This implies that the cascade described 



at the end of section 4.1 will also occur in the non-linear theory. 



What is different, however, is that the number of nodes in an eigenmode with given m is no 
longer fixed. Consider, for instance, the unstable modes of the 3-fold which come in a pair of 
opposite parity (see again Fig. Isj). The eigenmode has four nodes for Lp^ = ^ but eight for 



(fe = TT and 27r. A similar behavior is observed for the other n-folds (see Fig. [6]for the unstable 
modes of n = 4). This makes it difficult, in general, to assign the correct m to an eigenstate 
by inspection; its Fourier decomposition now displays a number of different frequencies (see 
Tab. [2] for the lowest modes of n = 2 and cpe = 2n). 

The natural resolution of this problem is to order the states with respect to increasing eigen- 
values. This provides a one-to-one correspondence which was also used to assemble Tab. [l] 
The zero modes, for example, are easily identified as the eigenvectors that go together with the 



vanishing eigenvalues. Our numerical calculations are consistent with the results of section 4.2 



one obtains the Fourier decomposition of three linearly independent eigenvectors each of which 
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m 


eigenvalue 


eigenmode 





1.41 


-0.2717 + 0.1075 cos (4s') - 0.0035 cos {8s') + 0.00008 cos (12s') 


1- 





0.1581 sin (s') + 0.3658 sin (3s') - 0.0127 sin (5s') - 0.0138 sin (7s') + 
0.0004 sin (9s') + 0.0004 sin (lis') 


1+ 





0.1581 cos (s') - 0.3658 cos (3s') - 0.0127cos (5s') - 0.0138 cos (7s') + 
0.0004 cos (9s') - 0.0004 cos (lis') 


2 





-0.3984 cos (2s') + 0.0199 cos (6s') - 0.0005 cos (10s') 


3- 


1.16 


0.3649 sin (s') - 0.1585 sin (3s') - 0.0292 sin (5s') + 0.0058 sin (7s') + 
0.0009 sin (9s') + 0.0001 sin (lis') 


3+ 


1.16 


0.3649 cos (s') + 0.1585 cos (3s') - 0.0292 cos (5s') - 0.0058 cos (7s') + 
0.0009 cos (9s') - 0.0001 cos (lis') 


4- 


6.08 


0.3988 sin (4s') - 0.0120 sin (8s') + 0.0003 sin (12s') 


4+ 


6.38 


-0.0760 - 0.3840 cos (4s') + 0.0115 cos (8s') - 0.0003 cos (12s') 


5- 


22.74 


0.0319 sin (s') + 0.00003 sin (3s') + 0.3975 sin (5s') - 0.00002 sin (7s') - 
0.0100 sin (9s') -0.0002 sin (lis') 


5+ 


22.74 


-0.0319 cos (s') + 0.00003 cos (3s') - 0.3975 cos (5s') - 0.00002 cos (7s') + 
0.0100 cos (9s') - 0.0002 cos (lis') 


6- 


56.67 


-44.4363 sn [1.8709s', -0.3003] + 43.7022 sin (2s') - 0.3299 sin (6s') + 
0.0033 sin (10s') 


6+ 


56.73 


0.0199 cos (2s') + 0.3984 cos (6s') - 0.0086 cos (10s') 



Table 2: Eigenmodes and -values corresponding to the modes m of the flat limit for n = 2 
and ipe = 2-71. In the table the scaled variable s' := ^^ = f is used. 
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can be written as a linear combination of the three zero modes $i,$j, and $fc- A fourth zero 
mode analogous to the one found in the small surplus limit does not exist. However, if isometry 
is relaxed in the numerical calculations, one finds a corresponding eigenmode which becomes 
unstable outside the small surplus limit. The Fourier series of this eigenmode displays exactly 



the same frequencies as the Fourier series (48) of the curvature k, i.e., m = {21 + l)n with 



/ G {0, 1, 2, . . .}. If the coefficient of the lowest frequency m = n \s set to the value of the 
corresponding coefficient of k, all other coefficients of the series of the mode are numerically 
found to be given by 2/ + 1 times the corresponding coefficient of the series of k. 

The eigenmode $o which corresponds to the constant mode (m = 0) in the small surplus 
limit now gets distorted by the Fourier mode m = 2n and small contributions of the multiples 
(4n, 6n, . . .). $o can be identified by looking at its eigenvalue: for (/^e — ^ the eigenvalue of 
$0 is given by n^ but diminishes if the surplus angle is increased (see again Tab. [l]). For the 
surplus angles under consideration the eigenvalue is always positive, the state (m = 0) is thus 
stable. 

One final observation concerns the eigenvalues: in the non-linear theory the two-fold degeneracy 
of certain pairs gets lifted. From Tab. [T]one identifies these pairs as the stable modes m = 4 
and 6 for the 2-fold, m = 6 for the 3-fold, and for higher values of m for the 4-fold. A quick 
glance in Tab. |2] provides an explanation in two distinct sources. One of these is the role of 
the constant term: whenever m is a multiple of 2n the eigenstate with even parity contains 
the constant term in its Fourier decomposition whereas the corresponding eigenstate with odd 
parity does not. The second is the isometry constraint: the Gram-Schmidt process breaks the 
symmetry between states of odd and even parity for those eigenstates with m = (21 + l)n. All 
remaining eigenstates come in pairs with equal eigenvalues, exactly as in the small surplus limit. 

5 Conclusions 

We have examined the stability of unstretchable sheets which are free to bend. The conical 
defect was studied in detail. Whereas cones with a surplus angle exhibit an infinite number of 
equilibrium states, all but the ground state with n = 2 are unstable. The dominant mode of 
instability will generally not be directly towards the ground state; instead there will be a cascade 
of instabilities. 

The fourth order differential operator, which arises in our analysis of stability, also possesses 
mathematical properties that are of interest in their own right. We have studied these properties 
to the extent that they touch on the issue of stability. 

Suppose that one removes a circular disk around the conical singularity at the center. The 
conical annulus obtained possesses additional isometric modes of deformation. These defor- 
mations can be decomposed into a conical part preserving the apex, as well as a part which 
destroys it. The resulting surface will be a more general tangent developable. Interestingly, 
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the linear operator determining the energy of deformations is not diagonal with respect to this 
decomposition. It is no longer even physically obvious that the two-fold conical ground state 
will continue to be stable with respect to isometric deformations destroying the cone. 

While we have limited our analysis to the behavior of an isolated conical defect, the framework 
developed here has more general validity. Analogous operators will occur in the study of the 
stability of related geometries: For example, for a d-cone subject to external constraints |3] the 



fourth order operator will have a form identical to Eqn. (34), with the potentials replaced by 
the relevant curvature function. Discontinuities at the points of contact, however, will make 
the analysis more delicate. The extension of our framework to sheets with a non-Euclidean 
metric [21, 22J is also possible but will be less straightforward: the analog of the operator 
C will be a partial differential operator; it will generally not be possible to integrate out the 
radial dependence. In these geometries one can exploit the rotational symmetry which predicts 
zero modes. When two or more defects are placed next to each other, however, one loses this 
symmetry and one would not expect these modes to show up. This makes the problem more 
challenging from a technical point of view but does not present an obstacle in principle. In any 
case, our analysis in this paper is a reliable point of departure for further research. 

We have benefitted from discussions with Martine Ben Amar, Osman Kahraman, and Norbert 
Stoop. We thank the Aspen Center for Physics as well as the Kavli Institute for Theoret- 
ical Physics for their hospitality. Partial support from DGAPA PAPIIT grant IN114510-3 is 
acknowledged. 

A First order deformations hKah and bK for an isometric sheet 

in general, the first order deformation of the extrinsic curvature is given by the covariant Hessian 
of the deformation vector projected onto the normal: 

bKa^ = -n-^a^^bX. (53) 

In an isometry, the trace K deforms by a divergence: 

6K = -n- V'^6X = -V„(n ■ V6X) . (54) 



The first identity uses Eq.(53) with Sg"''' = 0. The second involves pulling n inside the covariant 
derivative; the additional term generated, K'^^ea ■ Vb^X, vanishes for an isometry on account 



of the identity (12) 



It is convenient to decompose 5X into tangential and normal components: 

5X=V'"e„ + 0n. (55) 



Isometry (12) places constraints on these components 

VaV-fe + VfeV-a + '^Kab<P = . (56) 
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A simple calculation which makes use of the Gauss structure equations, W a^b = —Katn, then 
gives 

5ir = -V20 + V,(K'^V6)- (57) 

A useful alternative, non-divergence, expression for 5K is given by 

SK = (-V^ + 2Kg - K^)(p + {VaK) ^"^ . (58) 

Its counterpart for 5Kab can be expressed as 

5Kab = 5^Kab + 5\\Kab, (59) 

where 

5±Kab = -VaVb^ + KacK^cf), 
Sl^Kab = Kac^^r + KbcVar + r^cKab. (60) 



B Derivation of Eq.(29) 



Modulo two integrations by parts on the first term involving 5k" in Eq.([28j) we have 



-(fi:^)'$"-^[3/«2 + 2(l-q 



)]{k^)'^] 



We now use the local constraint (22) to eliminate the ^ dependence appearing in the last two 
terms {I, II) in favor of $. Note that the Euler-Lagrange equation £± = implies 

1, 



--[3K' + 2{i-q\)]K' 



Using Eqs.(22) and (61), it becomes possible to re-express the term // as 

-Lr„0 „/. ,^NT/ ^ ■. t ^ III ^ ^ -L 



2^2 



// = — [3fi;^ + 2(1 - q|)](K^)'<l> = k"^^ + -[3/?^ + 2(1 - C{)Y^ 



Modulo a total derivative. 
On the other hand 



fi:'"^$ = kk"$2 _ ^//^$/ ^ 



I = -(fi:^)'$" = k2$$" - k'^$" . 
We thus obtain (modulo total derivatives and the Euler-Lagrange equation) 

3 



/ + // = /t^$<l>"-^(/t'$')' + 



k" + -K-" + {1 - Ci\)k 



K<l>' 



^2$$" 



1 



1 



-KK 



k'Y + kN $ 



Summing terms, Eq.(29) follows 



(61) 

(62) 

(63) 
(64) 



(65) 
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